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Abstract: In this article the cross-section for the photodisintegration of a bound state 
is expressed, order by order in the multipole expansion, in terms of matrix elements be- 
tween states living on the three-dimensional torus. The motivation is to make the process 
amenable to Monte-Carlo simulations. The case of the deuteron is discussed. 
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1. Introduction 

A non-perturbative quantum field theory such as quantum chromodynamics (QCD) can 
presently be approached from first principles using Monte-Carlo methods, provided the 
theory is formulated in a finite spatial volume and in Euclidean time (see for example [1] 
for an introduction) . The former requirement is necessary in order to have a finite number of 
degrees of freedom to handle on a computer, and the second allows one to apply importance 
sampling methods. The spatial manifold is usually chosen to be a torus, which has the 
advantage of preserving translation invariance, at the cost of breaking the SO (3) rotational 
symmetry down to the discrete symmetry group of the cube. 

This setup is ideally suited to study the properties of the low-lying spectrum of the 
theory. However it is not immediately clear how to study scattering processes in this way, 
both because the scattering states are discrete in finite volume (with energy gaps that are 
in practice larger than the mass gap of the theory), and because the correlation functions 
require an analytic continuation before they can provide information on scattering phases 
and matrix elements involving timelike momenta [2]. In a numerical approach however, 
the analytic continuation represents an ill-posed problem, see e.g. [3]. 
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In a series of seminal papers [4, 5, 6], Liischer demonstrated that this limitation could 
be overcome for elastic scattering processes by relating the discrete two-particle spectrum 
on the torus to the scattering phases of the two particles (see also the older article [7], 
where the problem is addressed perturbatively in non-relativistic quantum mechanics). 
This analytic control over the two-particle spectrum allows one to extract other dynamical 
properties from stationary observables, in particular transition matrix elements from a one- 
particle state [8] or the QCD vacuum [9] to a two-particle scattering state. The derivation 
of the respective formulas given in these articles involves considering the maximal mixing 
between a one-particle state with a two-particle state on the torus under the influence of a 
perturbing Hamiltonian. However there are processes where a derivation along those lines 
is not possible (to the best of the author's knowledge), and the technical motivation of the 
present paper is to nonetheless derive a formula allowing one to extract a scattering matrix 
element from stationary observables on the torus in a different way. 

The process we look at is the photodisintegration of a bound state into a two-particle 
final state. Photoionization of an atom and photodisintegration of light nuclei are classic 
reactions in atomic and nuclear physics [10, 11, 12], and one reason we choose to study this 
type of reaction is because in these physics contexts it can be described in non-relativistic 
quantum mechanics, which is technically simpler. Indeed, except for the absorption of the 
photon, the particle numbers are conserved. We will generalize the quantum mechanical 
model (2.2) studied by Liischer on the torus to include the perturbative interaction of the 
particles with an electromagnetic field. The main technical novelty of this paper is the 
derivation of the relative normalization of a partial wave contributing to a finite-volume 
state and the same partial wave in infinite volume, Eq. (2.20). The main application follows 
from Eq. (4.4), which allows one in principle to calculate the photodisintegration of an s- 
wave bound state in the dipole approximation from numerical calculations on the torus. An 
important lesson one learns from the derivation is that the naive correspondence between 
the finite-volume states and the infinite-volume states must be examined individually for 
each process. 

Detmold and Savage investigated specifically the problem of calculating electroweak 
deuteron photodisintegration in lattice QCD several years ago [13]. They worked out a 
method in detail to determine the low-energy constants of the pionless effective field theory 
(EFT(^), [14]) relevant to magnetic dipole transitions. With these constants in hand, 
the photodisintegration process is predicted via the effective field theory. By contrast, we 
provide a method to directly compute the low-multipole matrix elements from lattice QCD, 
at the cost of having more stringent requirements on the volume of the torus. 

More broadly, the finite- volume methods to extract resonance properties have attracted 
a lot of attention in recent years, and it is appropriate to mention some of the work pub- 
lished recently. Ways to extend the formalism developed by Liischer to include the effects 
of multiple open channels have been proposed [15, 16, 17, 18]. In [19], interpolating opera- 
tors that are designed to couple to a resonance of given quantum numbers are analyzed in 
detail. At the same time there have been several recent numerical lattice QCD calculations 
of the interactions among pions in the p channel [20, 21, 22, 23], and of interactions in the 
two-baryon sector [32, 33]. 
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The outline of this paper is as follows. In section (2), we review the basic relations 
between the spectrum of states on the torus and the scattering phases and derive a relation 
between the wavefunctions in finite and infinite volume. Section (3) is devoted to coupling 
a two-particle quantum mechanics model on the torus perturbatively to photons, and the 
formula allowing one to calculate the photodisintegration process order by order in the 
multipole expansion is derived under stated assumptions. In section (4) the treatment is 
generalized so as to make contact with potential applications in QCD. We conclude on the 
prospects of such calculations and on future directions of interest. 

2. Relation between finite- and infinite-volume states 

Consider the discrete energy levels of the physical system consisting of two spinless particles 
of mass m on the L x L x L torus. We consider states in the A\ and in the T\ cubic 
representation. Under assumptions specified at the end of the paragraph, the energy levels 
are related to their infinite-volume scattering phase 5g(k) via [4, 5] 

kL 

6 l (k)+<f>(q)=nir, n € Z, 1= ^ ( 2A ) 

where (j)(q) is a known kinematic function 1 defined and investigated in Ref. [6], Appendix 
A. The lowest partial wave £ is and 1 respectively in the A± and T\ representations. The 
'effective momentum' k is defined from the energy of the state via the standard dispersion 
relation, E = |- in the non-relativistic case, fi being the reduced mass of the two particles. 
The relation (2.1) is strictly exact for a central potential that vanishes identically beyond 
a radius R (with R < L/2) and for angular momenta / > A, 

H = -±-A + Q A V(r). (2.2) 

For a wavefunction with spherical components Y>; m (r), 

oo £ 

V'W = E E Y im(0A)^im(r), (2.3) 
e=o m=-e 

the operator Q A acts as follows, 

a e 

Q A ^(r) = J2 Y em(0,(f>)Am(r). (2.4) 

£=0 m=-£ 

The relation (2.1) is valid for A strictly smaller than the second partial wave contributing 
in a given irreducible representation of the cubic group. Realistically, this means that it is 
a good approximation as long as the higher partial waves are small, which is guaranteed 
at sufficiently low energy. The relation generalizes in various ways. It remains correct 

x It is defined by tan 4>(q) = — ^J^i^ , where 2oo(b q 2 ) is the analytic continuation in s of Zao(s; q 2 ) = 
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up to terms exponentially small in L for a short-range potential (for example, a Yukawa 
potential). Remarkably, it also remains valid up to exponentially small terms in L in the 
relativistic case [5], except that k is then related to the energy level via the relativistic 
dispersion relation E = 2\/m 2 + k 2 . Finally, a generalization of (2.1) to arbitrary values 
of A exists [5]. 

It is interesting to consider how the energy level changes under a change in the box 
size. Eq. (2.1) implies 

Ldk(L) q<f/(q) 



k dL F E {k,L) 

where 



(2.5) 



F e (k,L) = k^^- + q( f> f ( q ). (2.6) 

The quantity Ff (k, L) that plays an important role in calculating the decay of an unstable 
particle on the torus [8, 9], appears repeatedly in the equations below. Eq. (2.5) provides 
a practical way of calculating it. 

On the other hand, one may consider a perturbation to the potential that leads to a 
change in the scattering phase A5i(k), for a given value of k. In a periodic box of fixed 
size, the allowed value of k must satisfy Eq. (2.1) both before and after the perturbation 
is added, for the same value of n. The difference of the two relations results in 

Ak 

A6f(k) = -F i (k,L) — . (2.7) 
This relation will be used in the next section. 

2.1 Relation between the finite-volume and the infinite volume wavefunction 

In this section we derive a relation between the finite-volume wavefunction \I/(r) in a 
given irreducible cubic representation 1Z and the wavefunction ip(r) of the infinite- volume 
scattering state, in the lowest partial wave i that contributes to 7Z. The starting point is the 
Hamiltonian (2.2) and we choose A such that Eq. (2.1) holds. Consider the wavefunction 
if}£ m in the lowest partial wave contributing to the given irreducible representation of the 
cubic group. In the interaction region r < R, there is a unique solution to the radial 
Schrodinger equation which remains finite at the origin and whose normalization is required 
to be such that r ip£ m = 1. This then determines the value and derivative of the radial 
wave function at r = R. Since there are exactly two linearly independent solutions to 
the radial Schrodinger equation, the matching conditions at r = R determine the two 
coefficients uniquely. Therefore the radial wavefunction in the lowest relevant partial wave 
is completely identical up to overall normalization to the corresponding infinite volume 
radial wavefunction. In the following we show how to calculate their relative normalization 
for a given normalization of the full wavefunctions. 

We denote a finite- volume wavefunction by \I/(r), and an infinite- volume wavefunction 
by ip(r). We assume that the former has unit norm, 

f d 3 r|^(r)| 2 = 1, (2.8) 
Jn L 
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where £Il refers to one periodic cell of dimensions L x L x L. In infinite volume, let the 
wavefunction take the form 

oo i 

V>(r) = E E Y em (9,<t>)^ m (r). (2.9) 
e=o m =-e 

Following the notation of [5], we write 

ip£m(r) = bi m ue(r; k), (2.10) 
with U((r; k) a solution to the radial Schrodinger equation 

I-" + k 2 ~ 2nV(r)) A m (r) = (2.11) 



normalized according to 

\imr~ e u e (r;k) = 1 (2.12) 
(in other words, bf m = lim r ^o r ^^m.( r ))- The general solution for r > R then reads 

u e (r; k) = a e (k)j e (kr) + f3 e (k)n e (kr) (2.13) 
for some coefficients ae(k) and (3g(k). The scattering phase bg is related to them by 



a e (k) + iMX) 
ae(k) - i/3e(k)' 

Now for each £ we choose the values of the bi m such that 



(2.14) 



b em (a t (k) + ift(fc)) = Smfn^/^e^ (2.15) 

for one particular value of m 6 { —£, . . . , +£} . The asymptotic r — > oo form of the radial 

wavefunction is then i^imir) = \j ^--gp sin(/cr — £^ — 5g). This implies that the states arc 
'energy- normalized', i.e. two states of energy E and E' in one particular partial wave £ 
normalized as in Eq. (2.15) satisfy (E = E' = 

/poo 
d 3 ri;(r;£,E)*^(r;£',E') = \b m \ 2 5 u , I dr r 2 u e (r; k)*u e (r; k') = S U 'S(E - E'). 

(2.16) 

If the interaction potential is varied (while remaining short-range), V — >■ V + AV, 
the discrete energy levels in the box will change. On the other hand, the phase shift also 
changes in infinite volume. These two changes have to be compatible with the Liischer 
relation (2.1) between finite- volume spectrum and phase shifts. 

The change in a non-degenerate energy level is given by first-order perturbation theory 
of quantum mechanics 2 , 

AE= [ d 3 r^(r)*Q A AV L (r)^(r). (2.17) 



2 On the torus, the relevant potential is Vl(v) = "Ylnez 3 ^ / (l r + n M) an< A similarly for AVl- We choose 
the angular momentum cutoff A between the lowest and the next-to-lowest partial wave in a cubic irrep. 
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On the other hand, the change in the phase shift is given by the generalized Born formula 
(see e.g. [24], parag. 133), which for our normalization of the wavefunction takes the form 

/"OO 

ASi=-7T r 2 drAV(r)\^ M (r)\ 2 . (2.18) 
J o 

The variation of E and of 5g must be related through Eq. (2.7) with Ak = AE and 
given by the dispersion relation. Combining (2.17, 2.18 and 2.7), one finds 

J™ r 2 dr A y(r)fe(r)| 2 = H jf dr *(r)*Q A AF L (r)$(r). (2.19) 

This equation is valid for any potential AV of range R < L/2. 

We will exploit in particular Eq. (2.19) in the T\ representation for I = 1 and A = 2. 
We can for instance choose a sufficiently short range potential AV, in which case one 
obtains the relation between the slopes at the origin of the two wavefunctions 3 . More 
generally we obtain in the I = 1 channel the relation (^e m (r) = j dQ Yi m (6, </>)* ^(r)) 

\4>iMr)\ 2 = H • ■ |*i,^(r)| 2 , (r < L/2). (2.20) 

To summarize, this equation confirms the statement that the finite- and infinite-volume 
wavefunctions are simply related to each other. The infinite volume, energy-normalized 
I = 1 wavefunction is proportional to the i = 1 component of the finite-volume, unit- 
normalized wavefunction in the T\ representation. The proportionality coefficient can 
be determined from purely spectroscopic measurements on the torus. Finally, we will 
choose the phase convention that V'io(^) and ^io(^) are both real. 

2.2 Higher partial waves 

In the region where the potential vanishes, the absolute normalization of the higher partial 
waves contributing to the finite-volume wavefunction can be determined. 

We consider again the T\ representation for £ = 1 and A = 2. In the 'outer' R < 
r < L/2 region, where V(r) vanishes, the wavefunction is the solution of the Helmholtz 
equation [5], 

V(r) = v 1 rnG 1 rn{r,k 2 ). (2.21) 
The radial function multiplying Y\fh(9, <j>) is [5] 

k 2 

*im(r) = -^«im [ni(kr) + Mim,im h(kr)} (2.22) 

On the other hand, we have seen that the radial wave function is proportional to the radial 
function in infinite volume. Comparing Eq. (2.22) with Eq. (2.20), one recovers Liischer's 
quantization condition, 

ai = Mlrn,lm{<l)Pl ( 2 - 23 ) 

3 Since the lowest partial wave dominates in any case at small r, the angular momentum cutoff is not 
needed in Eq. (2.19) if AV is sufficiently short-range. 
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as well as the condition that determines the normalization factor viq. Parametrizing gener- 
ically at = p^cosSt, jig = pgsmdi, one finds 



2ii d£4ir r . . 



For I > A the partial wave of the state ^(r) is given by 



®im(r) = --^vx^ Mifh/m je( kr ) = J F ^ k L ^ ksinSi Mifn,tm 3e(kr) (2.25) 
while the infinite- volume, energy-normalized wave function is given by 



2fik 

i>em(r) = 5 m rnd — je(kr) (2.26) 
since we neglect scattering in higher partial waves. 

3. A simple two-particle system and radiation 

We now couple the two-particle system discussed in the previous section to the photon. 
Much of this section is textbook material [10, 11], but serves as the preparation for the 
situation on the torus. Consider first, somewhat more generally, a non-relativistic iV-body 
system coupled to the electromagnetic field, 

N 1 

c=l c 

1 1 N 2 

= 2M ( P - I>cA c ) 2 + m Yl m ^ (Pc " tcA c ) -^- d ( Pd - e d A d )) ,(3.2) 

c,d=l 

where M = Y^=i m ci P = Yl^=iPc an( i we abbreviate A(r c ) = A c . We will restrict 
ourselves to the case of two particles, whose properties we index by a and b. Introducing 
as usual the center-of-mass and relative coordinates, 

m a r a + m b r b 

R= — , r = r b -r a (M = m a + m b ), (3.3) 

the Hamiltonian (3.1) can be written in the form 

tfkin = ^ (P - (e a A a + e b A b )f + i- (p - - ^A,))' (3.4) 

where tt = ™" a ™l is the reduced mass of the two particles. We note that the non- vanishing 
commutation relations are now [Pi,Rj] = \pi,Vj] = —iSij. 

We will be looking at one-photon processes and therefore only need the terms linear 
in the vector potential. We define 

Hl = ~ m {p > GaAa + ebAbh hl = ~\ {p ' { ^ Ab ~ ^ Aa)} (3 - 5) 
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where the curly braces represent anticommutators. Introducing a quantization box of size 
L q (which is unrelated to L and will be sent to infinity at the end), the gauge field is 
expanded in plane-wave eigenmodes, 



A(t, r) = J2 y ^ (a k ,,e a (k) e^~^ + a^fe)* e'^'^) , (3.6) 

with [afc j(7 , afc/ jCr /] = Skk'^ua' ■ For long wavelengths compared to the size of the bound state, 
we may expand the exponential. To calculate the anticommutator, we use the relation 
{pj, r{\ = ifJ-[H, TjTi] + eujLi, where L = r x p is the orbital angular momentum operator, 
and obtain 

e a ■ {p, e ik r } = 2p e a - fi[H, (e CT ■ r)(fe ■ r)] + iL ■ (k x e a ) + 0{k 2 ). (3.7) 

Thus 



^ = -^£^3 { aAk)e^ R -^\(% - ^^[H ,r ■ e„(k)} (3.8) 

-»\% + ^)lH ,(e a (k)-r)(k-r)] 

b a 

+ + %)L • (k x e a (k)) + 0(k 2 )) + h.c.}. 

In the same way, one finds 



H i = E \l^ q { a ° {k)etk ' R{e ° {k) • P) ( e « e_< " *" P + W^ k ' r ) + h-c-}- (3-9) 

In the center-of-mass frame, the matrix element of this term vanishes. 
3.1 Transition matrix element: dipole approximation 

We now want to study the cross-section for the photodissociation of a bound state of 
particles a and b. We assume for simplicity that this bound state is a scalar (pure s-wave 
bound state). 

Let us first assume that the wavelength of the photon is long compared to the size r s 
of the bound state, k^r s <C 1. Then the leading contribution to the matrix element of the 
interaction Hamiltonian in a power expansion in k^r s is 



(VfN^;fc 7 ,<T> = Uf -¥ x /^e ik ^ R (^--^-)[H ,r-e a (k 7 )} Vi) + 0(A; 7 ) (3.10) 



This expression determines the cross-section in the dipole approximation. Using a defi- 
nite angular momentum basis for the final state, we see that only its p-wave component 
contributes. According to Fermi's Golden Rule, the transition probability per unit time is 

dP 

— = 2nmh I \^)\ 2 p(E { ), (3.11) 
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where p(E) is the density of final states. For energy-normalized scattering states, this is 
just one. Calling 

J d 3 r Mr)*r^(r) = r fi , (3.12) 
the transition probability per unit time reads 

^ = ^(^-^)VMfer)T fl | 2 , (3.13) 

where we have used energy-conservation, Ef = E\ + w 7 . To obtain the cross-section, we 
must divide the transition rate by the photon flux. The latter is in the present case one 
per volume L 3 times the speed of light, and the cross-section for photodisintegration of the 
bound state is 

<r t =i = 4ttV(^ - ^) 2 u; 7 |e CT (fc 7 ) • r fi | 2 , (3.14) 
The key quantity to calculate is therefore the matrix element rg. 

3.2 Photodisintegration on the torus in the dipole approximation 

The treatment of the photodisintegration process laid out so far in this section carries over 
to the torus with little change. One difference is that the momenta accessible to the photons 
are discrete, fc 7 = ^n. The main question of interest here is whether the infinite-volume 
matrix element rg is accessible in the finite-volume theory. 

We work in the center-of-mass frame, P = 0, and assume the two particles to have equal 
masses. The initial bound state must therefore be moving with a momentum — fc 7 . In a 
non-relativistic treatment, the wavefunction describing the relative motion inside the bound 
state is independent of the total momentum, except for the torus boundary conditions. The 
spatial boundary condition of the internal wavefunction is either periodic or antiperiodic for 
equal-mass particles [25]. In the following paragraph, we focus on the internal wavefunction 
of the system. 

On the torus, both the initial and the final state internal wavefunctions contain an 
infinite number of partial waves. The former is in the A±, the latter in the T\ represen- 
tation 4 . However if the initial state is compact with a radius r s a few times smaller than 
the box size, then the only angular momentum component that is not exponentially sup- 
pressed in the volume is the s-wave component 5 [26, 27]. Since the position operator is 
a pure 1 = 1 operator, the only partial wave that can be reached is I = 1. Furthermore, 
since the position-space contributions to the matrix element are localized at r < r s , the 
matrix element rg would be the same as in infinite volume (up to exponential corrections), 

4 As far as the initial state is concerned, if the periodic/antiperiodic boundary conditions break the cubic 
symmetry, the wavefunction belongs to an irreducible representation of the reduced symmetry group. This 
does not affect the discussion. 

5 Indeed, the two linearly independent solutions to the free Schrodingcr equation arc and K t ,\. 

The latter however diverges at the origin, and must therefore be excluded for I > A. That leaves , 
which rises exponentially at large distances, say e KT . But since all the partial waves must be of the same 
order at r ~ L in order to satisfy the periodic boundary conditions, its coefficients must be of order e~ 2fiL , 
since the s-wave radial wavefunction is of order e~ Kh on the edge of the box. 
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if the normalization of the p-wave component of the final state were the same. With the 
normalization of states we have chosen, we find for the finite- volume matrix element Rr, 
using Eq. (2.20), 

Mfc 7 ) • r fi | 2 = i| • . K{ki) . i^p. (3 . 15) 

ah irk 

To repeat, the infinite-volume final state is meant to be a pure p-wave, while the finite- 
volume state belongs to the T\ representation. Eq. (2.20) - and therefore also Eq. (3.15) - 
assumes m a = mi,. Thus the cross-section for photodisintegration in infinite volume (3.14) 
can be expressed in terms of the finite-volume matrix element up to exponentially 
suppressed corrections. 

3.3 All-order multipole expansion 

If one does not expand in the momentum of the photon, the photodisintegration cross- 
section in infinite-volume is proportional to the square modulus of the matrix element 



hi 

I 



^■k 1 ,a)e- ik r R = J^-. (3.16) 
d 3 r Vf(r)*{ie CT (fc 7 ) • V, ^e^ fc ^ - ^e'^^ }^i(r). 



The standard multipole expansion of a vector plane wave with a wavevector in the z 
direction reads [28] 



(e x + iae y )e ikz = 



j e (kr)X^ a + cr^V x j e (kr)X £ja 



<r = ±l, (3.17) 



with C£ = i e y/4iv(2£ + 1). For the reader's convenience we summarize in appendix the main 
properties of vector spherical harmonics Yj^ (#,</>), of which the Xg m are a special case. 
Denoting /f (r) = -jJ= j e (kr)Y ej(T (8,<f>), one can rewrite 

{ex + i<TCl/ )c ite = - ^ ( " ctV (/f(»") + »" • - afe 2 r//(r) + irx V//(r)) = . 

e=i r Z£z 

(3.18) 

Now one can convince oneself that upon inserting the expansion (3.18) into the matrix 
element (3.16), each term labeled by the index i connects the initial bound state to a 
scattering state with orbital angular momentum given by L 2 \ip{) = £(£ + l)|^f). 

On the torus, the final state is chosen in the T\ representation. Its wavefunction 
contains infinitely many partial waves. Since we have related the normalization of the 
wave function in finite- volume to the normalization in infinite volume (see Eq. (2.20) and 
(2.25)), we can express the matrix element on the torus in the following way, 



(^—\ hl \^;k„a) = ^^^^^\ hl \^ tlhl a) (3.19) 

oo 

+ sin5i Y M ^M {^f ,m=rj) \ h Mii k ~ n °)), 
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Each term in this equation corresponds to a multipole amplitude. Furthermore, every term 
in the multipole expansion of the interaction Hamiltonian (3.16, 3.17) contributes exactly to 
one term in the series (3.19). From a practical point of view, the matrix elements calculable 
on the torus thus correspond to a linear combination of infinite- volume matrix elements to 
final states of definite angular momentum with calculable coefficients. These coefficients are 
L dependent, therefore by measuring the finite- volume matrix element at several L values 
and fitting them to a truncated version of formula (3.19), one can in principle determine 
multipole the matrix elements in several low-lying partial waves, provided such a truncation 
in £ is justified. 

We note that if the scattering phases for £ > 3 are not negligible, the coefficients 
of the series (3.19) change, but its structure remains the same, because in each partial 
wave, the radial wavefunction is proportional to the corresponding infinite-volume radial 
wavefunction. 



4. Potential applications in QCD 

In this section we investigate to what extent the results of section 2 and 3 can be applied 
to the two-nucleon system in QCD. 

4.1 Transition matrix element of a general electromagnetic current 

We now allow the charge and current density (p,j) to have a general form. The linear part 
of the interaction with the photon field then reads, in Coulomb gauge, 



hi 



+ H I ^H mt = - J d 3 xj(x)-A(x). (4.1) 



For j(x) = i ^{pa 5{x — r c )} one recovers the model of section 3. From that prepara- 
tory discussion it is clear that the relevant matrix element for the process of photodisinte- 
gration is 



J d^e^.^K^l^i). (4.2) 



The multipole expansion (3.17) of the photon plane wave can just as well be applied to 
this more general interaction. The lowest order term in ky comes from the gradient term 
in Eq. (3.18). Generically, when the photon field is a pure gauge, A = VG, the interaction 
Hamiltonian gives 

- J d 3 x j(x) ■ VG(aj) = -i H, J d 3 x p(x) G(x) . (4.3) 

In the present case G = ff{r) + r ■ V ff{r) is an eigenfunction of the operator L 2 with 
eigenvalue £(£ + 1). At low momentum transfer, the dominant term is £ = 1, and therefore 
the partial wave excited in is £ = 1. The relation between the finite- volume matrix 
element (4.2) and the infinite-volume matrix element at low energies is therefore the same 
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as in Eq. (3.15), 



(e=l,a) 



d 3 x e a (k 7 ) -j(x) e* 



(4.4) 



dA; Fx(k,L) 
~dE 



Tl,<7 



d 3 a; e CT (fc 7 ) ■ j(aj) e* 



+ 0(fc 7 ). 



We remind the reader that the normalization of the infinite-volume scattering state in this 
equation is {^\E)\^\E')) = 8 U > 6(E - E'). 

4.2 Magnetic transition near threshold 

We have so far assumed that the two particles are spinless and considered the transition 
from an s-wave bound state to a p-wave scattering state, which is the dominant one at 
low energies, because J = to J = transitions are forbidden according to the selection 
rules for one-photon processes [11]. However in the presence of spins, a transition t = to 
i = is possible, and in fact dominates at sufficiently low energies, as has been known for 
a long time. Fermi discovered in 1935 that the radiative neutron capture cross-section by 
hydrogen is non-negligible near threshold [29], and provided the correct explanation. The 
process is essentially a transition from a 2S+1 Lj = 1 Sq scattering state to a 3 Si bound 
state, the deuteron. The spins of the proton and neutron play a central role in this very 
low energy regime. In particular, the amplitude is non-vanishing only because the potential 
in the spin-triplet channel is different from the potential in the spin-singlet channel. This 
leads to s-wave radial wavefunctions of different energy that are not orthogonal to each 
other. 

The capture cross-section is related to the disintegration cross-section by Odis = fb-°"cap> 
where k is the relative momentum of the proton and neutron, and fc 7 is the momentum of 
the photon [11]. Here we will consider the disintegration process d + 7 — > p + n, which can 
be investigated on the torus along the same lines as above by including the spin degrees of 
freedom. We still ignore for simplicity the tensor force, which would lead in particular to 
a d-wave component in the deuteron. A contribution of the form 



le 
2" 



—cr a x [p a , 5(x - r a )} + — a b x \p b , S(x - r b )} 
m a m b 



3m{x) 

to the electromagnetic current, with the corresponding interaction Hamiltonian 



J dx j M {x) ■ A(x) = ~[ —cr a ■ B a + —a h ■ B t 



(4.5) 



(4.6) 



2 \ra a m b 

induces (at lowest order in the multipole expansion) a magnetic dipole transition. Here fx a 
and jjL b correspond to the magnetic moments of the proton and neutron in Bohr magnetons, 
e.g. 2.78 for the proton. Denoting the two-nucleon spin states by x(S, Tn s ) (for instance, 
X(!, !) = xtxt and X(0,0) = -^(XaX b ~ XaX b )) the initial state with m = 1 is of the 
form ip = ^i(r)x(l, 1) and the final state ip f = V'f ( r )x(0> 0)- Decomposing the current into 
an isovector (<r n — a p ) and an isosinglet (<r n + a p ) part, one easily checks that only the 
former makes a contribution to the 3 Si — > 1 Sq transition 6 . 

6 Yet twisted boundary conditions cannot straightforwardly be used to calculate this cross-section at very 
low momenta, for the reasons given in [30]. 
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As in the case of the 3 <Si — > °Pi transition driven by the electric dipole operator, 
the 3 Si — > 1 Sq matrix element on the torus can be related to the infinite- volume matrix 
element, 



(£=0) 



I 



d 3 x e CT (fe 7 ) ■ j(x) e 1 ' 



dk F (k,L) 



d 3 x e a (k r )-j(x) e l 



(4.7) 



dE irk 

4.3 Probing the nucleon-nucleon potential 

We return to relation (2.17), which relates the energy shift due to a change in the inter- 
particle potential to the wave function for a non-degenerate state. If we choose again the 
angular momentum cutoff A between the lowest and next-to-lowest partial wave allowed 
by the cubic symmetry and choose the perturbing potential to have a range smaller than 
L/2, the relation is equivalent to, up to exponentially small corrections, 



(4.8) 



where £ is the lowest partial wave contributing to the given cubic representation (I = and 
1 respectively in the Ai and T\ representations), and m £ {— £, ■■■+£}. A determination 
of AE thus gives access to information on the infinite-volume wavefunction. For instance, 
choosing the one-parameter family of Yukawa potentials, 



AV(r) = g 2 



-fir 



4nr 



(4.9) 



the Laplace transform of rtp£ m (r) 2 is obtained by varying fj,. An inversion of the Laplace 
transform allows one to evaluate the unperturbed potential via 



(pe(r) r z 



(4.10) 



We note that the absolute normalization of the wavefunction is not required to determine 
V(r). In the case of two nucleons, let us consider then a channel with a definite total 
spin S and definite orbital angular momentum I. We disregard here the effects of the 
tensor force which mixes several orbital angular momenta into a state of definite total 
angular momentum J, leaving this issue to future investigation. It may be unrealistic 
to invert the Laplace transform numerically, but the wavefunctions obtained by solving 
phenomenological nuclear potentials can be inserted into the right-hand side of Eq. (4.8) 
and compared to the left-hand side calculated in lattice QCD Monte-Carlo simulations. In 
this way, Eq. (4.8) provides a consistency check between lattice QCD and phenomenological 
approaches to the two-nucleon system. We expect this basic idea to remain valid when the 
local potential V(r) is replaced by a non-local one. 

There is no simple way of implementing an angular momentum cutoff on the perturbing 
potential. In general the energy shift (4.8) then corresponds to summing over £ the term 
appearing on the right-hand side. In order to have a simple interpretation, the method 
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must be applied in the regime where the effective momentum k is small in comparison with 
the range of both the internucleon potential and the perturbing potential AV(r). In this 
way the contribution of higher partial waves is kinematically suppressed. 

How can one implement the method at the quark level? The quarks can be given a 
Yukawa coupling to a massive, color-neutral scalar boson. The nuclcon mass itself will 
then shift by order g 2 , and this effect must be subtracted in order to consider the change in 
binding energy of the two-nucleon state in Eq. (4.8). An attractive feature of this method 
to probe the internucleon potential is that only spectroscopic 'measurements' are involved . 
However, very accurate numerical data is required. 

5. Conclusion 

We have derived a relation, Eq. (4.4), which provides a recipe to calculate numerically the 
amplitude for the photodisintegration of a scalar bound state. The result is based for one on 
the relation between the lowest partial wave contributing to the finite- volume wavefunction 
and the infinite-volume wavefunction in the same partial wave; they are identical [5] up to 
overall normalization, and the latter is given by Eq. (2.20). And secondly, the result relies 
on the multipole expansion of the transition amplitude. 

In order to comment on the prospects of exploiting relation (4.4) for the photodisinte- 
gration of the deuteron in lattice QCD, we recall the assumptions made. An important one 
is that the initial bound state is well contained in the box. This assumption ensures for 
instance that the matrix element of the dipole operator between the initial bound state and 
the final scattering state is entirely due to the p-wave component of the latter, and that 
this contribution is not affected by the 'cutoff' in the radial variable at r = L/2. Secondly, 
we assumed that only the p-wave phase shift is significant, which is only justified at low 
energies. This assumption can be relaxed thanks to the general analysis performed in [5]. 
Thirdly, the box size must be large enough for the smallest momenta of the photon, i.e. 
the momentum transfer between the initial and the final bound state, to be in the relevant 
low-energy region. 

The requirements on the physical volume of the torus are thus very strong, since the 
deuteron is a very weakly bound state, -Ebind ~ 2.2MeV and the energy range where the 
dipole approximation works well is below lOOMeV. A box size of lOfm or more is required. 
Approaching the problem from heavier quark masses, where the scattering phases have 
more 'natural' sizes, could be possible if the deuteron remains bound, as suggested by re- 
cent lattice calculations [32, 33] (some earlier effective field theory estimates [34, 35] had 
suggested that the deuteron unbinds quite soon as the light quark masses are increased 
from their physical value). Finally, apart from the kinematic requirements, obtaining an 
acceptable signal-to-noise ratio in the two-nucleon channel is a tremendous challenge in lat- 
tice QCD, see for instance [36, 32, 33]. As an exploratory calculation, it may be interesting 
to apply the methods proposed here in (quenched) lattice QCD at heavy quark masses 
and very large volumes as in [32] . Perhaps the process can also be investigated in effective 

7 This is to be contrasted with a recently proposed procedure to determine off-shell nucleon potentials 
that depend on the interpolating operator, see [31]. 
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theories along the lines of [37, 38]. Whether there is a choice of boundary conditions that 
alleviates the strong requirements on the box size deserves further investigation. In any 
case, we hope that the techniques presented in this paper will be useful to address further 
scattering processes. 

Studying the disintegration of nuclear bound states through operators other than the 
vector current is also of great interest. As already emphasized in [13] in the lattice QCD 
context, the axial current is a particularly important case, since it is relevant to the disin- 
tegration of the deuteron induced by neutrinos (see for instance [39]). This reaction was 
instrumental in determining the flavor of the neutrinos arriving from the Sun in the SNO 
experiment [40]. 
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A. Vector spherical harmonics 

For the reader's convenience we briefly describe the properties of the vector spherical har- 
monics and how the multipole expansion Eq. (3.17) comes about, following the treatment 
of [11]. The vector spherical harmonics Yff^O,^) are vector functions of the (9,4>) angles 
which are simultaneous eigenfunctions of the operators J z and J 2 , 

J z = L z + S z , L = -irx V, Si = e { x (.) (A.l) 

and 

J 2 = Jl + J 2 + J 2 . (A.2) 

The eigenvalues are respectively M and J(J + 1). In addition, the Yjf x have parity (— 1) . 
Since L and S commute, they are constructed by the standard addition law of angular 
momentum, 

l l 

*jfl(M)= E J2 C a (J,M;m,m')Y em (9,<j>)e m ,, (A.3) 

m=— £ m'=— 1 

where 

ei = -={e x + ie y ), e = e z , e_i = —= (e x - ie y ). (A.4) 
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The vectors X^ m appearing in Eq. (3.17) correspond to the special case 

X lm {0A) = Yui(0A) = -j=== LY tm {9,<f>). (A.5) 

The vector spherical harmonics are orthonormal, 

J dn YjM(9, cP)* Ytf^O, 4>) = 5jj,5 MM ,S ee , (A.6) 

Since the vector spherical harmonics form a complete set, an arbitrary vector field can be 
expanded as 

oo I 

Mr) = E E [ftm(r)Y^ + gtm(r)Y% +1A + h em {r)Y%_ ltl ] (A.7) 

£=0 m=-£ 

The radial functions fi m ,gtmi he m are obtained from the integrals 

f tm [r) = J <mY%(n)*-A(r), (A.8) 

9tm(r) = J dQ Y e m e+hl (n)* • A(r), (A.9) 

h tm (r) = J dQ Y%_ ltl {ny • A(r). (A.10) 

The commutation relation [Vx, J] = implies that the curl of an expression with 
some definite value of J and M must be a linear combination of vector spherical harmonics 
with the same J and M. In particular, for an arbitrary radial function f(r), we have 

V x (f(r)X em ) = f + (r)Y% +1A - f-(r)Y$_ 1A , (A.ll) 

where f±(r) are related to f(r). Thus the expansion (3.17) of the plane wave is of the 
form (A.7). Since only m = ±1 appears on the right-hand side of Eq. (3.17), each term 
is individually orthogonal to e z . Finally, the specific radial coefficients f,g, h of Eq. (3.6) 
can be obtained from the integrals (A.8-A.10). 
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